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Abstract — For fixed c Prolate Spheroidal Wave Functions i\) n ,c form a basis with remarkable prop- 
erties for the space of band-limited functions with bandwith c and have been largely studied and 
used after the seminal work of Slepian. Recently, they have been used for the approximation of 
functions of the Sobolev space H s ([— 1, 1]). The choice of c is then a central issue, which we address. 
Such functions may be seen as the restriction to [—1,1] of almost time-limited and band-limited 
functions, for which PSWFs expansions are still well adapted. To be able to give bounds for the 
speed of convergence one needs uniform estimates in n and c. To progress in this direction, we push 
forward the WKB method and find uniform approximation of ip n ,c in terms of the Bessel function 
Jo while only pointwise asymptotic approximation was known up to now. Many uniform estimates 
can be deduced from this analysis. Finally, we provide the reader with numerical examples that 
illustrate in particular the problem of the choice of c. 
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1 Introduction 

Traditionally, the prolate spheroidal wave functions (PSWFs) have been used for solving various 
problems from physics and signal processing, see [7J [TU [HJ [22 . Perhaps, the PSWFs were specially 
known for their important contributions in solving many problems from antenna theory, see for 
example [14]. Note that for a given real number c > 0, called bandwidth, the PSWFs denoted 
by {tpn c('))™>0j were known as the eigenfunctions of the Sturm-Liouville operator's L c defined on 
C 2 ([-l,l])by 



LM = (1 - x*)pt - 2x^ - c 2 x^. (1) 
ax ax 



The eigenvalues (— Xn(c))n>o ar e fixed by the requirement that the eigenfunctions ipn,c(x) are 
bounded as \x\ — > 1~ . To the best of our knowledge, in [15] . C. Niven was the first, in 1880, to 
give a remarkably detailed theoretical, as well as computational study of the eigenfunctions and 
the eigenvalues of the above Sturm-Liouville problem. In their pioneer work [T21 EH 1221 1221 Hi] > 
D. Slepian, H. Landau and H. Pollak have shown various important properties of the PSWFs and 
their associated spectrum. Among these properties, they have proved that the PSWFs are also the 
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eigenfunctions of the compact integral operators F c and Q c , denned on L 2 ([— 1, 1]) by 

FcW>)(*) = [ Sm f X ~? My)dy, Qc{f)(x) = [ e lc *yf(y)dy. (2) 

J -i n \ x ~ y) J-i 

As a result, they have shown that the PSWFs exhibits the unique properties to form an orthogonal 
basis of L 2 ([— 1, 1]), an orthonormal system of L 2 (R) and an orthonormal basis of B c , the Paley- 

Wiener space of c— band-limited functions defined by B c = |/ g L 2 (R), Support / C [— c, c]| . 
Here, / denotes the Fourier transform of /. The PSWFs are normalized by using the following rule, 

/ \i/> n , c (x)\ 2 dx = l, ( \ifj n . c (x)\ 2 dx = n>0. (3) 

J-l JM An(C) 

(A rt (c)) rl is the infinite sequence of the eigenvalues of F c , arranged in the decreasing order 1 > 
A (c) > Ai(c) > ••• > A n (c) > ••• . 

Recently, there has been a growing interest in the study of the quality of the spectral approx- 
imations by the PSWFs and the building of PSWFs based numerical schemes for solving various 
problems from numerical analysis, see [3J |31 [5J [55]. In particular, in [3], the author has shown 
that a PSWFS approximation based method outperforms in terms of spatial resolution and stability 
of timestep, the classical approximation methods based on Legendre or Tchebyshev polynomials. 
The authors of [5] were among the first to study the quality of the approximation by the PSWFs 
in the Sobolev space H s ([—1, 1]), s > 0. In particular, they have given an estimate of the decay of 
the PSWFs expansion coefficients of a function / 6 H s ([—1, 1]), see also [4]. Recently, in [25], the 
author studied the speed of convergence of the expansion of such a function in a basis of PSWFs. We 
should mention that the methods used in the previous three references are heavily based on the use 
of the properties of the PSWFs as eigenfunctions of the differential operator L c , given by JT]). They 
pose the problem of the best choice of the value of the band- width c > 0, for approximating well a 
given / g H s ([—1, 1]), but their answer is mainly experimental. It has been numerically checked in 
[H [25] that the smaller the value of s, the larger the value of c should be. 

When starting this work, our aim was to clarify the way c should be chosen. We also wanted 
to use the properties of the PSWFs related with the Fourier transform. These ones are not easy 
to use for numerical computation because of the rapid decrease of the eigenvalues of Q c , but they 
are nevertheless fondamcntal. To choose both c and the index of the partial sum used in the 
approximation, the main difficulty is that this requires uniform bounds for the PSWFs for c variable. 
The first main contribution of this work is to derive a whole set of qualitative and quantitative 
properties of the PSWFs. Some of these results were already known in the literature from mainly 
numerical evidences and/or asymptotic and heuristic justifications. 

More precisely, the asymptotic behaviors of eigenvalues A„(c) and functions ipn,c are well-known 
in different conditions, for c tending to or oo for instance. These kinds of behaviors are confirmed 
by numerical evidences [501 123] ■ We give rigorous proofs, which are valid for the whole range of 
values of n and c, or, for most of them, under the condition that c 2 < Xn(c), which is asymptotically 

C 7T 1 

equivalent to the condition given in [4], that — < 1, where c* = — -(n + — ). An emblematic estimate 

c n ^ 2 

that is frequently used in the literature is the following one, 

|</Vc(l)| < + (4) 

This is stated in [3D] without a rigorous proof, but justified by asymptotic expansions and numerical 
evidence. This is a key point to prove that the sequence A„(c) decreases very rapidly. An analytic 
proof of estimate ([4]) for all n and c seems out of reach. Nevertheless, by pushing forward different 
methods, we obtain uniform estimates on the values of ip n ,c that are sufficient to obtain an exponen- 
tial decay of A„(c). We have a gain compared to these authors in the sense that we have a complete 
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proof, while we loose on constants that appear in the exponential decay compared to theirs, which 
are comforted by numerical evidence. The two approaches are clearly complementary, and we do 
hope that having rigorous proofs valid in the whole range of c and n adds to the understanding of 
the expansions in PSWFs. 

Our second main contribution is to use different properties of the PSWFs as eigenfunctions of the 
differential operator L c and the finite Fourier integral operator Q c and give a study of the quality of 
approximation by the PSWFs in the Sobolev spaces H s ([— 1, 1]). More importantly, our study tries 
to give a satisfactory answer to the previous important problem of the choice of the parameter c. To 
this end, we are led to study various inequalities and asymptotic results associated with the PSWFs. 

This work is organized as follows. In Section 2, we prove various uniform estimates of the 
PSWFs. In particular, we use the WKB method and provide an explicit uniform approximation 
over [0, 1] of the PSWFs. Also, we derive several useful bounds of the values the PSWFs and their 
first derivatives. Moreover, we give two results on the approximation of the PSWFs by Legendre 
polynomials. In Section 3, we give some practical and useful estimates of the decay of the eigenvalues 
and the Legendre coefficients associated with the PSWFs. Our results improve the existence decay 
results given in (6] [26]. In Section 4, we first give the quality of approximation by the PSWFs in 
the set of almost time and band-limited functions. Then, we combine these results with those of 
Section 2 and give a new L 2 ([— 1, 1])— error bound of approximating a function / e H s ([— 1, 1]) by 
its Nth. terms truncated PSWFs series expansion. The proof of this bound is based on the use of 
the quality of approximation of almost bandlimitcd functions by the PSWFs. Also, by using the 
Fourier characterization of periodic Sobolev spaces, we give another approach for functions that 
extend into periodic functions with the same regularity. These new estimates provide us with a 
way for the choice of the appropriate bandwidth c > to be used by a PSWFs based method for 
the approximation in a given Sobolev space H s ([— 1,1]). In Section 5, we first give an overview 
concerning the computational aspects of the PSWFs, then we provide the reader with two methods 
for the computation of the PSWFs approximate expansion coefficients of a function / G H s ([— 1, 1]). 
Finally, we give some numerical examples that illustrate the different result of this work. 

2 Uniform estimates for the PSWFs 

In this section, we prove that the PSWFs ip n ^ c are uniformly approximated in terms of the Bessel 
function Jo when c is not too large, namely, when c 2 < aXn( c )> where a < 1. The proof of this 
approximation result is based on the use of the well known WKB method for the study of asymptotic 
behavior of the solutions of the perturbed differential equations of the Sturm-Liouville type, see for 
instance [T7j. We emphasize the fact that the asymptotic behavior that we shall describe is well- 
known for a single value of the variable and n tending to infinity (see for instance [U [5] and the 
references there). What is new here is the possibility to push forward the methods in order to 
have uniformity. This allows us to have uniform bounds for the functions ip n c within the prescribed 
constraints on the parameters and recover partially the bounds given in [5D] from numerical evidence. 

2.1 Uniform approximation of the PSWFs by the WKB method 

2 

In the section we note q n (c) — x c / c \ ■ Most of the results of this section are obtained under the 
relatively weak constraint that q n (c) < 1. Also we always assume that "0n,c(l) > 0. We prove 
here that ipn,c can be approximated by a simple explicit function (once one knows Xn( c ))i U P to a 
multiplicative constant A(n, c) that takes into account the normalization that we have chosen for 
■i/Vi.c- This constant is not explicitly known, but tends to 1 when q n (c) tends to 0. On another side, 
we compute explicitly the dependance in q of constants and do not allow the notations O and o, 
except for functions of one variable and when there is no ambiguity. 
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To simplify notation, we now skip the parameters and note tp, q and \ n . We only leave the 
parameter n for the eigenvalue Xn so that the asymptotic behavior of each term is clear. With these 
notations, -0 satisfies the differential equation 



d 



— [(1 - X 2 W(X)] + Xn(l - qX 2 )^{x) = 0, xe [-1, 1]. 



(5) 



Because of the parity of the PSWFs, we can restrict to the interval [0, 1], which we will do from now 
on in this sub-section. We then note 



S(x) := S q (x) = 



l-qt 2 



(It 



(G) 



which defines a homeomorphism on the whole interval and is smooth in [0, +1). It is well-known as 
an elliptic integral of the second kind. We look for -0 under the form 



i/)(x) = cp(x)U(S(x)), <p(x) = (1 - .x 2 )- 1/4 (l - qx 2 y 1/4 . 
Replacing S' by its explicit value, Equation ([5} becomes 
d 



d;T' 



-(1 - x 2 ) 1/2 {l - qx 2 ) 1/2 LpU' + (1 - x 2 )tp'U + Xn(l - qx 2 )ifiU = 0. 



(7) 



(8) 



One can easily check that terms in U' disappear with this choice of ip, and the equation satisfied by 
U on the interval [0, +1) may be written as 



with 

Let us write 
so that 



U" + (Xn + hi)U = 0, 



h!{S(x)) := ^-'(l - qx 2 )' 1 -^ [(1 - **W(x)] . 



(9) 



Q(x) := (l-x 2 )(l~qx 2 ), 



(10) 



< P '/<P = - l Q'/Q. 

It follows that hi o S is a rational function with poles in ±1 and ±J^, which may be written 



h 1 oS=—(l-qx 2 y 1 
Id 



Q! 



(1-^)1^) -4-((i-^)^ 



Q! 



Only the pole 1 is of interest (because of the assumption on q) . One easily computes the residue at 
1, which is a — —4(1 — q)^ 1 ■ So h\(S(x)) can be written as 



hi(S(x)) 



(l-g)" : 
8(1 -x) 



h 2 (S(x)), 



with /12 o S a rational function without poles on [0, 1]. Elementary explicit computation proves that 
this function is bounded by 2(1 — q)~ 3 . 

Next we want to replace (1 — x)^ 1 by a function of S(x). This possibility will be a consequence 
of the following lemma. 
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Lemma 1. One has the following inequalities, valid on the interval [0, 1]. 

2(l-«)(l-a0 <S 2 q (x)<4(l~x). (11) 

At one has 

1 < S,(0) < |. 

Moreover extends into a holomorphic function in a neighborhood of 1 and takes the value 

2(1 — q) at 1. Finally 



< S q (x) - y/2(l - g)(l -*)<\ |* ^ 1/2 ■ (12) 
Proof. The first inequality comes from the facts that the function is decreasing and J^ 1 -yj= = 

. ,2 

2yl — x. Moreover the square root of ~_jj f extends into a holomorphic function in a neighborhood 
of 1. Extending it into an entire series in 1 — t, multiplying it by y/1 — t and taking the antiderivative, 
we find that 

1—x 

elementary bound 



3/2 



we find that S ^ also extends into an entire series in 1 — t. For the last inequality, we use the 



1 + t V 2 - (l-g)V2- 



By dividing the previous inequalities by yl — t and then integrating the different sides of the in- 
equalities over [x, 1) and using ©, the estimate p^|) follows at once. □ 

It follows from Lemma [T] that 2 (i-q) T^x ~ Six)' 1 ex te n ds into a holomorphic function in a neigh- 
borhood of 1. Moreover it is bounded by rrj^p • ^° we can W1 "it e this difference as h^(S(x)), with 
/13 a bounded and continuous function on [0, 5(0)]. If we define F := — (hi + fta/4), we have proved 
that F is continuous and satisfies 

< (j^a- (13) 

This is summarized in the following lemma. 

Lemma 2. For g < 1 there exists a function F :— F q that is continuous on [0, 5(0)] and satisfies 
(1131) such that U is a solution of the equation 

U"{s) + (xn + ^) U(s) = F(s)U(s), s e [0, 5(0)]. (14) 
Before using the properties of such an equation, let us translate for U the fact that ip has norm 

1. 

Lemma 3. We have the equality 

1 

(l-ga; 2 )- 1 |J7(5(x))| 2 d5(x) = 1. (15) 



In particular, 

,S(o) 

(1-q) < 2 / \U(s)\ 2 ds < 1. (16) 
./o 

Proof. Just use the expression of -0 in terms of U and the expression of 5'. □ 
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Let us come back to Equation (fT4|) . which we will consider as an equation in U with second 
member FU. It is well known that the associated homogeneous equation has the two independent 
solutions 

Ui(a) = xl/^MVx^s), U 2 (s) = xTV^Y (^s), 

where Jo (resp. Yq) denotes the Bessel function of the first (resp. second) type. The Wronskian of 
Ui and U% is given by 

W{U u U 2 ){s) = Ui^U^s) - U[(s)U 2 (s) 

= sXn [Jo{Vx^ s ) Y o(Vx~^ s ) - J'o(Vx7i s ) Y o(Vx~^s)} 
2^ 



when using the well-known identity for the Wronskian of Bessel functions. Then the general solution 
of HH), according to the method of variations of constants, leads to the following identity. 



U(s) 



AU^s) + BU 2 (s) + 



x / V^[Jo(Vx7 t s)Y {y^t) ~ J (^t)Y (^s)]F{t)U(t)dt. 



(17) 



Let us pose 



K n (s,t) := VstXn [Jo(\ / Wis)Ya(^/x^t) - Jo(Vx^t)Y (y/x^s)] 



Using the fact that 



supVs(|Jo(s)| + |*o(s)|) < 00, 

s>0 



we find that the kernel K n is uniformly bounded. Using Schwarz Inequality, f|16[) and the estimate 
on F, we find that 



- / K n (s,t)F(t)U(t)dt 



< 



Cs 1 ' 2 



(18) 



Here C is a uniform constant, that does not depend on q. 

In particular, this integral remains bounded when divided by s 1 / 2 for s tending to 0. Let us 
prove now that B = 0. Since tp extends continuously to 1, the function 



U(S(x)) 



tr = iJ{x) x ((1 + x){1 9x2)1/4 



remains bounded for x tending to 1 and the same is valid for U(s) / y/s for s tending to 0. On another 
side, Jo(s) remains bounded while Yo(s) is not bounded for s tending to 0. This forces B to be 0. 
So (|17l) can be rewritten as 

U{s) = A{ y fx^s) l ' 2 J Q { y / X ^s) + -^— f K n (s,t)F(t)U(t)dt. (19) 

The following lemma is needed in the proof of the main result of this section, given by Theorem 



Lemma 4. Let I = 

and c, such that 



Xn Jo 



t (Jo(t)) 2 dt. Then there exists a constant C independent of n 



25,(0) 



< 



a 



(20) 
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Proof. We first recall that 1 < 5,(0) < tt/2. It is well known (see [2]) that 



t(Jo(t)) 2 d< = y [(J (x)) 2 + (Ji(x)) 2 ] , x > 0. 



Moreover, we recall the asymptotic behavior of Bessel functions, see [T§] 



sup | v / a;<A' (ir) | < c vi SU P 

x>0 a;>0 



J "( x ) - \l — cos ( x - 2 " ~ 4 



<^2' ">-!• (21) 



Here, c v ,d v are constants depending only on the order v. We easily get 



That is 



Ji(VXnS 9 (0)) 



1" v/X^^gC ) 



< 



A' 



Xn(^(0)) 2 



A = doco + dici + 2 



25 9 (0) 



<-S=, C = 2(d c + d 1 c 1 )+4J-. 



(22) 
□ 



We are now able to state our main theorem. We mention explicitly the parameters. In the 
statement, we assume that ^ n)C is such that VVi.c(l) > 0. 

Theorem 1. There exist constants C, C with the following properties. Assume that the parameters 
n, c are such that q = c 2 /Xnip) < 1- Then one can find a constant A :— A(n, c) < M such that, for 
< x < 1, 

M*) - ^ n(c l?^y^f' (a;)) + RnM (23) 



(l-.T 2 ) 1 /4(l_ g:E 2)l/4 



with 



C 



SUp \RnA x )\ < CgXn(c) 1/2 , C q = ^yj. 

xe[o,i] (l-g) 1 ' 5 /" 



(24) 



Moreover, for any integer n satisfying Xn > (^ 25 (0) 
constant A is bounded above and below as follows, 



, where C is as given by Lemma 4> the 



Sk. 

v'TT 



25 9 (0) 



Proof. Using (|19J) . we can write 



7T C" 

2S,(0) ^ 



< A < 



1 + V2- 



Xn 



25 9 (0) h __2r Si 



(25) 



with 
where 



R n , c (x) = (1 - a; 2 )- 1 /^! - qx 2 )- 1 '^^)), 



R(s) 



2\/Xn(c) A 



K n {S(x),t)F{t)U(t)dt. 
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To bound R n (x), we proceed as follows. From (TT81 . we have sup ^ < ^ — Also, 

*e[o,i] y/Six) (l-'7) 3 v / X^ 

from (ED, we get ^sup ^ (1 _ x2)1/4(1 _ g , 2)1/4 < (^^J- M ° reOVCT ' SmCe 

i?„, c (x) = - , 2 )-V4 (1 - 

then from the previous discussion, we conclude the bound of |i? n , c (x)|. To prove (|25p . we first pose 
V(s) = ( v / X^s) 1/2 ^o(VX^s) for s G [0, 5,(0)]. Note that E7 = AV + R. We deduce from the previous 
estimates that the L 2 norm of R is bounded by v / 2C 9 x„(c) -1 / 2 . Let us recall the inequalities (|16l) . 
It follows that 

v/T^g- V2C, X n(c)- 1/2 < |A|||V|| 2 < l + V2C qXn (c)- 1/2 . (26) 

To conclude it is sufficient to prove that the norm of V is conveniently bounded above and below. 
After a change of variable, we have to use the bounds above and below for the quantity /, given by 



Lemma 4. More precisely, by writing A 2 ||V||2 = A 2 



25,(0). , 25 9 (0) 



-) + 



and by using (|20|) and 



flU, one gets (J25J). □ 

Theorem Q] may be used on the whole interval [—1, +1], using the parity of the functions ipn,c- 
It gives precise indications on the functions ip n . c , both quantitatively and qualitatively. We could 

as well have written the bounds for the rests in terms of n^ 1 instead of Xn(c) -1 ' 2 - Indeed, from the 

inequality n(n + 1) < x„(c) < n(n + 1) + c 2 one deduces that 

/ \ n(n + 1) , 
n(n + l)<Xn(c)< \_ ■ ( 27 ) 

The next corollary gives bounds for tp n ^ c . It is an immediate consequence of Theorem [I] and 
the fact that | Jo{s)\ + s 1 / 2 \ Jq{s)\ is uniformly bounded for s > 0. We use the same notations as in 
Theorem [TJ 

Corollary 1. There is a constant C such that, for C q = C(l — <z)~ 4 ; the two following inequalities 
hold. 

sup N <! < Cgxlic (28) 

sup^^l - z 2 ) 1 / 4 ^^)! <C q . (29) 
Remark that the first bound is sharp since 

^„ >c (l)=A(n,c)xn(c) 1/4 . (30) 

Next we state as a lemma the fact that, for q close from 0, the constant A is close from 1. 

Lemma 5. Let a < 1 and < K < 1. Let C be as defined by Lemma 4- There are constants 

(ji 

H\ = H\{a,K) and iJ 2 = H2(ot,K) such that, for q < a and n satisfying < K\/l — a, the 

constant A(n,c) in Theorem^ satisfies the inequality 

\A 2 {n,c)~l\<H 1 {a,K)q + H 2 (a,K)xn(c)- 1/2 . (31) 
As a consequence, under the same assumptions on q, 

ip n , c (l) 2 -n- ~ <H 3 qn + H 2 . (32) 



Proof. We first note that 

*4^<a, ( o)<f 

By using the previous inequalities as well as (PZ51) . we get 

(1 - g) 1 | gT^p < A» < (1 - g) -V 2 ^ c , fl ffi/. - ( 33 ) 
We find in particular that 



l-A 2 < q +(C' + 2^ 2 C a ) X - 1/2 



a 



Under the condition that —— < Kyi — a, which implies that the denominator on the right hand 
side is larger than 1 — K, we also have 

A 2 - 1 < + (1 - - a)" 1 / 2 (2 3 / 2 C a + 2^) X - 1/2 . 

V 1 — a 

By combining the two estimates, one obtains (|3"TT) with iJi = (1 — a) -1 / 2 and i?2 = (1 — K)^ X {1 — 

a) -l/2 (2 3/2 CQ+2C 2 ) _ 

Next, we use (f5T|) and the inequality 



Xn(c) 1/2 -n-i < 



2 " ^Xn 
to conclude that 

1 



< p 2 -i|vx^+IVx^-(« + 1 /2)l 

□ 

Remark that (|52")1 can be compared to the conjecture that the left hand side is negative. Even 
if much weaker, it will nevertheless be sufficient to prove the exponential decrease of the sequence 
A„(c). 

2.2 First approximation by Legendre polynomials 

We now discuss the possibility to modify x«( c ) m Formula (1231) . It is easily seen that Theorem Q] 
remains valid when x„(c) is replaced by Xn(c) + /?, for some fixed constant. Indeed, one just has to 
change F into F + f3. This is particularly interesting for c = 0, where, instead of Xn(0) one can take 



(n + 1/2) 2 . With this modification, the well-known fact that P n (l) = y n + | leads to the fact that 
^4 = 1. Also, 5(x) = arccos(x). So, we have the following corollary of Theorem [TJ 

Theorem 2. The normalized Legendre polynomials satisfy the following, uniformly for < x < 1. 

1/2 



= \jn + \ ( T"" 8 ^ ) Jo((n+l/2)arccos( a; )) + 0^) (34) 

The uniformity of the approximation seems new, even if the formula is not. 

We can also replace x«( c ) by (n + 1/2) 2 in Formula ([23)) for all c, except that now the constant 
C has to be multiplied by (1 + c 2 ). If we want that the second term behaves like a rest, we have to 
assume that c 2 is sufficiently small compared to Xn(c) x / 2 . This leads us to the condition: 
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• Condition £(a, e): 

So let us assume this condition and replace Xn(c) by (n + 1/2) 2 in Formula (|2"5)) . It is easy to see 
that the rest of the proof remains unchanged. We do not give too many details because we will use 
another method later on. So there there exists some constant B = B(c,n) such that we can write 



= B (1 _ x 2)i/4 (1 _ ga; 2)i/4 + (35) 

with |i? n . c | < Cx n (c)~i~ e , with C a uniform constant depending only on a and e. Moreover (we 
also proceed as in Lemma [3]), the constant B is such that 

\B{n,c)-\\ <H X n(c)-*- e . 

So we can replace it by 1, which generates an error of the order Xn(c)~^~ e ■ In order to go further, 
we want to replace S q (x)) by arccos(cc) without too much loss, in order to be able to prove that 
ipn,c is close from the Legendre polynomial. Since S q {x) > (1 — q) 1 ^ 2 arccos(ir), the difference 
arccos(x) — S q (x))is also bounded by Hx n (c) ~^~ e arccos(a;). This implies that | Jo((n+^) arccos(x)) — 
Jo(( n + ^)S q {x))\ is bounded by Cx„(c) _e . The other factors are also easily shown to be close to 
the analogous terms for the Legendre polynomial. We have proved finally the following proposition. 

Proposition 1. Under the condition C(a,e), there exists C such that 

sup \ip n , c (x) -P~ n {x)\ < Cn~ 2t . 

Different authors have proposed a second term as a candidate for the expansion of ip n ,c in terms 
of Legendre polynomials for q tending to zero, see for instance Boyd's paper [4]. Our proof cannot 
lead to this because our rest is not small compared to this second term. 

The previous proposition does not answer a natural question, which is the approximation by P n 
for c tending to 0. We will give another method for this. 

2.3 Second approximation by the Legendre polynomials 

Proposition 2. Let a < 1. There exists a constant M a with the following property. For all n and 
c > such that q n {c) < a, then 

i — i c 2 

sup Vn.cfo) - P n (x) \< M a =, (36) 

xe[-i,i] 01 + 1/2 

Proof. The proof follows the same lines as the proof of Theorem[TJ We now see the equation satisfied 

by ip := i)n,c as 

(1 - x 2 )ip" - 2xip' + n(n + l)Vw = (n(n + 1) - Xn + c 2 x 2 ) tp. 
The homogeneous equation 

(1 - x 2 )^" - 2xtp' + n(n + l)tp n , c = 

has as solutions the Legendre polynomial P n on one side, the Legendre function of the second kind 
Q n on another side. Moreover (see [1]) the Wronskian is given by 

W(P n ,Q n )(x) = Y-^. 



10 



The function G :— n{n + 1) — Xn + c 2 a; 2 is uniformly bounded in modulus by c 2 . By the method of 
variation of constants we can write 



ip{x) = AP n (x) + BQ n (x) + 



n + 1 /2 



(1 - y 2 ) (P n (x)Q n (y) - P n (y)Qn(x)) G(vW{y) dt. 



We pose 

L n (x, y) = (1 - y 2 ) (K(x)Qn(y) - ^(y)0»(as)) • 

We claim that L n is uniformly bounded, independently of n, x, y. Let us take this for granted and 
go on for the proof. Then, for some uniform constant H, the error term satisfies 



1 



n + 1/2 



L n (x,y)G(y)ip(y)dy 



< H 



n + 1/2' 



The consideration of the behavior of each term when x tends to 1 implies that B is equal to and 
we have, for again H a uniform constant, 

_ r 2 

\4>-AP n \<H 



n+1/2 ' 



Since both functions if) an( A Pn have L 2 norm 1, we have the inequality 



\l-A\<H 



n + 1/2' 



Consequently, we have 



sup - P n (x)\ < sup \ip n , c (x)-AP n (x)\+\A-l\ sup \P n (x)\ 

xe[o,i] xe[o,i] xe[o,i] 

r 2 

< M n 



□ 



It remains to prove that L n is uniformly bounded. It is a consequence of the fact that yl — x 2 \ P n (x) \ 
and Vl — x 2 \Q n (x)\ are bounded independently of x and n. For P n , it is a consequence of (|29|) . For 
Q n it is given by the following lemma. 

Lemma 6. There exists a constant H such that, for all n, we have the inequality 



sup y/l - x 2 \Q n (x)\ < H. 
xe[-i,i] 



Proof. We recall that Q n (x) = -P„(x)log 



1 + a: 
1-x 



Vn-l(x) = - 



P n {x)-P n {t) 



x — t 



dt = ^ n + 1/2 



— V„-i(x), where 

1 Pm-l{t) P n -m(t) 



E 



\ m ^Jm - 1/2 y/n - m+ 1/2 



We can restrict to x > 0. We recall that, by (|29|) . (1 — x 2 ) 1 ^ 4 \P n (x)\ is uniformly bounded. So the 
required inequality follows at once for the first term. Next we consider V n -i. Using again (1^1) . it is 
sufficient to prove that 



E 



y/n + 1/2 



i x my/m— l/2 x /n — m+ 1/2 

is uniformly bounded, which is elementary: just consider separately the cases m < n/2 and m > n/2. 
This finishes the proof of the proposition. □ 
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2.4 Uniform bound estimates of the PSWFs and their derivatives. 

In this paragraph, we give various estimates of the bounds of the i/Vi.c an( i tp'n c over [ — lj !]• Moreover, 
we give bounds for the successive derivatives of the PSWFs at x = 0. The content of this paragraph 
will play later on a major role for the approximation of polynomials by their series in PSWFs. By 
using the results of Theorem 1 and its corollary, we prove the following proposition that provides us 
with a first set of uniform bounds. 

c 2 

Proposition 3. There exists a constant C depending only on a < 1 such that, for q = — < a, 

Xn 

then 

sup W{x)\<C Xn 5/ \ (37) 
xe[-i,i] 

sup (1 - x 2 )\i>'{x)\ <C^[. (38) 
xe[-i,i] 

Proof. To prove ([37]) . it suffices to consider the case x E [0, 1]. Then by using the identity 

(1 - x 2 )^'(x) = Xn ! (1 - Qt 2 )m dt, (39) 

J X 

one gets 

(1-X 2 )\yj'(x)\ <Xn(l-X 2 ) [ \i,(t)\dt. 

J x 

Finally, inequality (f28|) gives us (|37]l . 

To prove (J35J), we start again from ([23) ■ We let s = S(t), t e [0, 1]. Here, S(t) is given by @. 
Let L(s) be the function defined on [0, 5 g (0)] by 

//i _ *2U/2\ 1 / 2 

L(s)= L (s-H)=(i- q t 2 )^ ^ 1 s y> j . 

It is easy to see that L(s) is of class C 1 on [0, S q (Q)], with a norm bounded by a constant that 
depends only on a. Moreover, by using ([2"3"]l . we get 

l r s(x) ,1 

(1 - qt 2 )4>{t) dt = AxH A / L{s)sJ (^s)ds+ R n At)dt. (40) 



Since (sJi(s))' = sJq(s), then an integration by parts gives us 



S(x) 



XT j L(s)sM^s)ds = [L(5(x)) v ^) x y 4 v ^)Jo( v ^r5(x)) 

1 



'Xn 



Six) 

L'(s) Xn /i sJ 1 (^s)ds. (41) 



\JXn Jo 

By combining $5§§, (gDJl, dH} and using (|ST]t. we easily conclude for ([55]). □ 

Note that the constants in the previous proposition degenerate for q tending to 1. We will see 
below than one has nevertheless uniform bounds in a neighborhood of for \ip n ,c\ and \ij)' n c \. This 
is the subject of the following proposition. 

Proposition 4. There exists a constant C such that, for all n > and c > 0, 

sup {l-x 2 ) l '^ n , c {x)\ < (2 X n(c)) 1/4 . (42) 

B6[-l,l] 

sup (l-x 2 M' n Jx)\ < C(c 2 +xn(c)) 3/4 . (43) 

x£[-l,l] 
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Proof. We shall again use the differential equation ([5]) satisfied by ipn,c, which we note only ip in 
this proof, and the identity 



+ 1 



(1 - X 2 M( X )\ 2 dx + C 2 / X 2 \^{ X )\ 2 dx = Xn(c) 



(44) 



The above identity is obtained by multiplying ([5]) from both sides by i/j(x) and then integrating the 
result over [—1, 1]. We can consider only x > and replace 1 — x 2 by 1 — x. By a simple integration 

by parts applied to the quantity 2 / (1 — t)ip(t)ip'(t)dt, one obtains the following identity 

J X 



|V>(t)| 2 * + 2 / (l-t)ip{t)^'{t)dt. 



(45) 



We prove that the left hand side of the previous identity is bounded by (1 — x) 1 ! 2 y^Xn(c)- To this 
end, it suffices to bound the second integral. This is done as follows. 



(l-t)Tjj(t)i(>'(t)dt 



< 2vT 



< 2yT 



li>{t)ip' (t)dt 
1 i 1 . - r r' 

(i-t 2 W(t)\ 2 dt 



m)?dt 



1/2 



(46) 



The last inequality is a consequence of (|4"4"1) . By combining (l4"5"j) and (T4"6")l . we conclude for (|42p . 
Next, we write 



(l-xY\^'{x)Y = -2 / {l-t)\i)'{t)\ 2 dt + 2 (l-t)VW'(*)* 



(47) 



Again, we use ([5]), multiply this later by \J\ — t 2 and get 



(i - t 2 f 2 r(t) = 2W1 + (ch 2 - Xn )Vi -t 2 m, t e [-1, 1]. 



Since \t\, Vl — t 2 < 1 and by using Minkowski inequality, one gets 



That is 



||(1 _ t 2 f' 2 r\\2 < 211^1^ *V||2 + (C 2 + Xn)||^|| 2 . 

U (l-t 2 ) 3 W(t)\ 2 dtj <2VxW)+Xn(c)+c 2 . 



(48) 



Again, by using Schwarz inequality in the second integral of (|47l) and then combining (j47]) and (|48l) , 
one gets g2]). □ 

The following proposition gives us interesting bounds for the successive derivatives of the PSWFs 
at x = 0. 

Proposition 5. For am/ integers n, fc > 0, satisfying k(k + 1) < Xni w £ /save 



1^(0) < (v^) K 1^ 



n.c v w y I 7 



for n even and k even, and 



mm < (^ fe_1 K,c(o)i , 



(49) 
(50) 
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for n odd and k odd. In particular, under the assumption that q — — < 1, there exists a constant 

Xn 

C, depending only on q and such that for any positive integer k satisfying k(k + 1) < Xn, we have 



<>(0) <C(^) k . 



(51) 



Proof. We first study the case where n = 2m, k = 2i are even integers. We show that for a fixed n, 
V>J$(0) has alternating signs, that is ^c(°) = (-l) fe/2 V>H(0) or equivalently ^(O)^"^ (0) < 0. 

Note that by an iterative use of ([5]), one can easily check that the ip^l(0) — ip^ k '(0) are given by the 
following recurrence relation, 



^(*+2)( ) = (k(k + 1) - x„)V (fe) (o) + k(k - i)c 2 v (fe_2) (o), k > 0, 



(52) 



with t/>(0) > 0, V (2) (0) = -Xn^(0). Note that V> (2) (0)^(0) < 0. By induction, assume that this is the 
case for the order k - 2, that is (0)V> (fe-2) (0) < 0. Multiplying ([52]) from both sides by ^ (fc) (0), 
using the assumption that k(k + 1) < Xn as well as the induction hypothesis, one concludes that the 
induction assumption holds for the order k. Consequently, we have, 



Let 7 fc 



^(*+3)(q) =(x„-ife(*+l)) ^ (fc) (0) +k{k-\)c 2 V (fe_2) (0) 
^( fc )(0) = m k (^/x^) k . Then jS3]) is rewritten as follows, 
k(k+l) 



k>0. 



rrik+2 



m k + k(k-l) — m fe _ 2 , 1-0(0)1 = mo- 

Xn 



(53) 



(54) 



From ?M 2 )(0) = -Xn^(0), we have to 2 — mo. Moreover, from (154"]) . we get 



fc(fc + l)\ ,„ a 

mk+2 < ( 1 I TOfc + fc(fc + 1) TOfc-2- 

Xn / Xn 



A simple induction gives us 

m k+2 < (l - Mi±!2(i _ g)J m < m = |^(0)|. 

This concludes for (|49p . Similarly, one can easily prove (1501) . Finally, by combining 
(Unj), cone concludes for (|ST|) . 



and 
□ 



3 Decay estimates for the eigenvalues and the Legendre co- 
efficients 

3.1 Decay estimates for the eigenvalues 

We first use the previous estimates to prove the exponential decay of the sequence A„(c). More 
precisely we prove the following. 

Theorem 3. Let S > 0. There exists N and K such that, for all c > and n > max(iV, kc), 

A„(c) < e - S{n - Kc) . 



— . It is a corollary of the following technical proposition. 



This can be compared with the numerical evidence that one has super-exponential decay from 
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Proposition 6. Let a < 1. Then there exists constants Mi, Mi such that 



\n(c) < M^C-f" 



c 2 c 2 
— exp(M 2 ^) 



(55) 



Proof. We start from the well-known equality, see [201 124) . that for any positive integer n, we have 
A„(c) = A' x A", with 

c 2n+1 (nl) 4 



A' : 
A" : 



2((2n)!) 2 (I> + 3/2)) 2 

(Vy r (i)) 2 -(n + l/2) 



exp 2 



(h 



(56) 
(57) 



Let us hrst consider A". We will use (1321) and ([36 



It follows that 



\(^ T (l)?-{n + 1/2)| < M{l + r 2 /n) 

< Mt 2 

X" <e 2M (e 2M ^Y . 



r> 1, 

T < 1. 



On the other hand, we have 



2((2n)!) 2 
ri/n 



< I - 



n i+£ 



(58) 



(59) 



Since - Y^ log(l + j/n) > Y/ log(l + a;) dx = 2 log2 - 1, then 



< e -(41°g2-2)n y?1 > L 



By using the previous inequality together with (|58|). one gets the required inequality ([55]) . □ 

Remark 1. Numerical evidence, see f20f . indicates that (ipn,T) 2 — (n + 1/2) < 0,Vt > 0. If we 

c ( ec \ 2n 

accept this assertion, then we observe that the sequence A n (c) decays faster than — y~J so that 
the exponential decay has started at [ec/4]. 

3.2 Decay estimate of the Legendre expansion coefficients 

In this paragraph, we study some decay estimates of the Legendre coefficients Recall that 

ipn,c{x) — YJ f3%Pk(x), V x € [—1, 1], with /3fc = / Pk(x)ip n , tC (x) dx. It is well known, see [Theorem 
fe>o 



3.4, [26]], that for a fixed positive integer n, we have 

2 1 



Vfe > 2([ec] + 1). 



(60) 



lM4c)|2 fc: 

Here fJ, n ( c ) is the eigenvalue of the operator Q c defined in the introduction, for the eigenfunction 
Tpn.c- It is related to A„(c) through the identity 

A «(c) = ^-|Mn(c)| 2 . 
The following lemma improves the decay estimate given by (|60[). 
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Lemma 7. For any real c > and any integer n > 0, we have 



Pk{x)lp n ,c{x) dx 



< - 



2) T(fc + 3/2) | M „(c)| 



, Vfc>0. 



(61) 



Proof. We start from the following identity relating Bessel functions of the first type and Legendre 
polynomials, see [2] 

11 e ix vp n {y)dy = i n X [ 1 ^J n+k {x), V x ? 0. (62) 



Since the ip n ,c are eigenfunctions of Q c , that is 



e lcxy ip n , c {y) dy = fi n (c)ip n<c (x), 



(63) 



then, by combining f)63[) and (|62[) and by using Plancherel Theorem, one gets 



1 



ft = / P k (x)^ n , c (x)dx 

-1 Mn(c) 
1 ' X 



1 /.l 



e lcxy P k {x)dx\^ c {y)dy 



Mn(c) 



vfe + i/2w — Jk+i/2(cy)i> n ,c(y) d v- 

V cy 



By using the previous equality together with the well known bound of the Bessel function, see [5] 



\J a {x)\ < 



2«r(a + l) 



, Va>-l/2, Vie 



(64) 



one gets 



< 



< 



y/ir(2k + l)c k 
2*+i/2r(fc + 3/2)| A t n (c)| y_! 



\y k UnAv)\dy 



^Ti{2k + l)c k 



2 fe + 1 /2r(fc + 3/2)|/i„(c)| V 2fc + l 



\?PnAy)\ 2 d y 



1/2 



1 



2 J r(fc + 3/2) |/x„(c)| ' 



□ 



Remark 2. Our decay bound given by &61]) outperforms in two directions, the one given by A60\) . 
Firstly, it is valid for any integer k > 0. Secondly, our decay bound decays much faster than the one 
given by i60\) . In fact, Since T(s + 1) > \/2~7rs s+1 / 2 e~ s for all s > 0, see f]}/, then it is easy to see 
that 



(!) 



1 



< 



1 



T(k + 3/2) |/i n (c)| - V2 fc+1 V 2fc + I J V 2 ^ 1 lMn(c 



1 1 



,Vfc>2([ec] + l). 



4 Quality of the spectral approximations by the PSWFs 

In this section, we first study the quality of approximation of almost band-limited functions by the 
classical PSWFs ij} nc that are concentrated on [—6,6], for some b > 0. Then, we extend this study 
to the case of periodic and non periodic Sobolev space H s ({— 1, 1]), s > 0. 
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4.1 Approximation of almost time and band-limited functions 

In this paragraph, || • || 2 denotes the norm in L 2 (K). We show that the set {ip n ,c( x ), n > 0} is well 
adapted for the representation of almost time-limited and almost band-limited functions, which are 
defined as follows. 

Definition 1. Let T = [— a, +a] and Q = [—6, +6] be two intervals. A function f, which we assume 
to be normalized in such a way that \\fW2 = 1, is said to be et— concentrated in T and en— band 
concentrated in Q. if 

\f(t)\ 2 dt<e 2 T , ±- j |/H| 2 ^<4. 

Up to a re-scaling of the function /, we can always assume that T = [—1,1] and ft = [— c, +c], 
with c := ab. Indeed, for / that is ex— concentrated in T = [—a, +a) and en— band concentrated in 
= [—6, +6], the normalized function g(t) — y/af(at) is ex— concentrated in [—1, +1] and cq— band 
concentrated in [—ab,+ab]. 

Before stating the theorem, let us give some notations. For / an L 2 function on K, we consider 
its expansion / = X^n>o a n'4'n,c in L 2 {[— 1, +1]), with, due to the normalization of the functions ipn,c 
given by ([3]), the following equality holds, 

[ +1 \f(t)\ 2 dt = J2K\ 2 . (65) 

We call Sjv.c/) the iV-th partial sum, defined by 

SN,cf(t) ■= 2J a„VVc(£)- (66) 

n<JV 

We write more simply Sat/ when there is no ambiguity. In the next lemma, we prove that Sjyf 
tends to / rapidly when / belongs to the space of band-limited functions. This statement is both 
very simple and classical, see for instance [2TJ[22] or Theorem 3.1 in [25] . 



Lemma 8. Let f € B c be an L 2 normalized function. Then 

c + l 



\f-S N f\ 2 dt<X N (c). (67) 

Proof. Since the set of functions -0„ jC is also an orthogonal basis of B c , the function / may be written 
on R as / = J2 n >o a n^n,c, with 



j |/(t)| 2 dt = ^|A n (c)|- 1 |a n | 



n>0 

The two expansions coincide on [—1, +1], and, from (1551) applied to / — Sjyf, it follows that 

r + l 



f |/ - S N f\ 2 dt < sup |A„(c)| V lA^c)!- 1 



a n \ 2 - 



i>N 

We use the fact that the sequence |A„(c)| decreases and to conclude. □ 
Next we define the time-limiting operator Pt and the band-limiting operator IIo by: 

PtU){x) = xr(x)f(x), IIn(/)(aO = [ e ixu f(cj)dw. 

The following proposition provides us with the quality of approximation of almost time- and band- 
limited functions by the PSWFs. 
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Proposition 7. If f is an L 2 normalized function that is ex— concentrated in T = [— 1, +1] and 
en — band concentrated in CI = [— c, +c], then for any positive integer N, we have 

/ f+l X 1 / 2 , 

(/ |/-5 JV /| 2 dtJ < en + VAiv(c) (69) 

and, as a consequence, 



\\f-FrS N fh<€ T + en + y/XN{c). (70) 

More generally, if f is an L 2 normalized function that is ct— concentrated in T — [— a, +a] and 
en — band concentrated in Q — [—b, +b] then, for c = ab and for any positive integer N, we have 



11/ - PrS N , c , a fh < ct + en + VAjv(c) (71) 
where Sjv,c,a flwes i/ie AT-t/i partial sum for the orthonormal basis -y=ijj n c (t/a) on [— a, +a]. 

Proof. We first prove (ffi9")) by writing / as the sum of IIq/ and g. Remark first that \g — 
SNg\ 2 dt < \\g\\2 < en- We then use LemmaOfor the band limited function Hnf to conclude. The 
rest of the proof follows at once. □ 

Remark 3. Let f be a normalized L 2 function that vanishes outside [— 1, +1] and is in H S (M). 
Then f gives an example of 0- concentrated in [— 1,+1] and e c -band concentrated in [— c, +c], with 
e c <M f /c s . 

Indeed, take 
which is finite by assumption. 

4.2 Approximation by the PSWFs in Sobolev spaces 

In this paragraph, we study the quality of approximation by the PSWFs in the Sobolev space 
H s ([— 1, 1]). We provide an L 2 ([— 1, l])-error bound of the approximation of a function / £ H s ([—l, 1]) 
by the N— th partial sum of its expansion in the basis of PSWFs. 

To simplify notation we will write / := [—1, 1]. In this paragraph, |) • H2 denotes the norm in 
L 2 {I). We should mention that different spectral approximation by the PSWFs in H s = H S (I) have 
been already given in [H |6j [25] . More precisely, the following result has been proved in [6] . Here 
a k{f) ■= J_ 1 f(x)ip k ,c(x)dx. 

Theorem 4. (Theorem 3.1 in 0). Let f € H S (I), s > 0. Then 

M/)l < c 6v- 2/3s ||/||h 

where C, 5 are independent of /, N and c. 

In |25j , the author has used a different approach for the study of the spectral approximation by 
the PSWFs. More precisely, by considering the weighted Sobolev space H r (I), associated with the 
differential operator 

V c u = -(1 - x 2 )u" + 2xu' + c 2 x 2 u, 




18 



and given by 

H r (I) = I / € L 2 (I), ||/|||„ = \\V r J 2 f\\ 2 = ^(XkYlfkl 2 < +™ 

[ fe>0 

the following result has been given in [25] . 

Theorem 5. (Theorem 3.3 in [25 ). For any / € H r (I), with r > 0, we /iai>e 

||/-5 W /|| 2 < {XN{c))- r ^\\f\\ Rr <N- r \\f\\ Sr . 

It is important to mention that the error bounds of the spectral approximations given by the 
previous two theorems, do not indicate how to choose a convenient value of the bandwidth c to be 
used to approximate a given / <E H S (I). By a simultaneous use of the properties of the PSWFs as 
eigenfunctions of the differential operator L c and the integral operator Q c , we give a first answer to 
this question. This is the subject of the following theorem. 

Theorem 6. Let c > be a positive real number. Assume that f £ H S (I), for some positive real 
number s > 0. Then for any integer N > 1, we have 

||/ - S N f\\ 2 < K(l + c 2 )- s / 2 ||/|| H . + Ky/XMcj\\f\\ 2 . (72) 

Here, the constant K depends only on s. Moreover it can be taken equal to 1 when f belongs to the 
space Hq(I). 

Proof. To prove (|72p , we first use the fact that for any real number s > 0, there exists a linear and 
continuous extension operator E : H S (I) -> H S (R). Moreover, if / £ H S (I) and F = E(f) £ H S (R), 
then there exists a constant K > such that 

\\F\\ LHm < K\\f\\ 2 , ||F|| h .(r) < K\\f\\ H s. (73) 

We recall that the Sobolev norm of a function F on M is given by 

\\F\\h m :=^j^(l + |^| a )*|/(0l a de- 
In particular, for F c— bandlimited, one has 

iiF|i| 2(R) <(i+ c 2 r s ii^ii^( R) - 

Next, if T denotes the Fourier transform operator and if 

Q = F-\F ■ 1[_ C)C] ), n = F-\F ■ (1 - l[_ CiC] )), 

then Q is c— bandlimited and F = Q + H. Moreover, since ||5||l2( R ) < ||F|| L 2( R ) and ||K||_j,2( K ) < 
c ~ s ||-P 1 ||ff=(R): tnen by using ([75]). one gets 

\\G\\ L 2 m < K\\f\\ 2 , \\H\\ 2 <K(l + c 2 r s / 2 \\f\\ H s. (74) 

Finally, by using the previous inequalities and the fact that Q is c— bandlimited, one concludes that 

\\f-s N f\\ 2 < \\g-s N g\\ 2 + \\n-s N Hh 
< V^N(c)\\g\\ L 2 m + \\H\\ 2 



< ^/)^K\\f\\ 2 + K(l + c 2 )- s \\f\\ Hs . 

This concludes the proof for general /. When / is in the subspace Hq(I), one can take as 
extension operator the extension by outside [—1,1], so that the constant K can be replaced by 
1. □ 
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Remark 4. This should be compared with the results of [251, given by Theorem 5. This has the 
advantage to give an error term for all values of c, while the first term in (1721) is only small for c 
large enough. On another side, Wang compares this specific Sobolev space with the classical one and 
finds that 

n/ii^<c(i+c 2 r /2 ii/ii^. 

For large values of N we clearly have xn(1 + c 2 ) <C (1 + c 2 ) -1 , but it goes the other way around 
when xn and 1 + c 2 are comparable. So it may be useful to have both kinds of estimates in mind for 
numerical purpose and for the choice of the value of c. 

Remark 5. Remark that the second error term in \12^ is negligible as soon as N is comparable to 
c, due to the super- exponential decay of the sequence Xn{c), and also to the fact that the L 2 norm 
may be very small compared to the norm in H s . A convenient choice of the value of the truncation 
order N c to be used is given by 

N c = min{N G N, V%v(c)||/|| a < OI/IIW- 

If N < N c , then the error bound in {72ty is concentrated in the quantity yj Ajv(c)||/||2 which can 
be significantly reduced by considering larger value of N. On the other hand, if N > N c , then the 
error bound is concentrated in the quantity c~ B \\f\\H s - Hence, larger values of N will not reduce in 
a significant manner this error bound. 

It may be useful to consider in particular the subspace H^ er of functions in H s (I) that extend 
into 2— periodic functions of the same regularity. For such functions, one can also use the norm 

ii/ii^ er .=E( 1+ (M 2 ri6fe(/))i 2 - 

Here, 

h(f) = ^= £ f(x)e-™ k *dx = j=f(kn) 

is the coefficient of the Fourier series expansion of /. A precise error analysis of the quality of 
approximation by the PSWFs in the space Hp er is given by the following theorem. We should 
mention that the proof of this result is essentially based on combining the Fourier characterization 
of Hp er with some properties of the PSWFs as the eigenfunctions of the finite Fourier transform 
operator Q c . 

Theorem 7. Let s > 0, c > 0, be any positive real numbers and let f G Hp er ([— 1, 1])- Then for any 
integer N > 1, we have 



\\f-S N f\\ 2 < (1/2+-) £ ||^c|| 2 A„(c)||/|| 2 + C - s ||/-/ [c/7r] ||^ cr . (75) 
V C n>iV 

Here, f[ c /n] * s the truncated Fourier series expansion of f to the order [— 1 . In particular, for any 
positive integer N satisfying q — c 2 /\N < lj we have 

11/ - S N f\\ 2 < K q Ij2 VA^A„(c)||/|| 2 + c- s \\f - / [c/w] ||H. ar> (76) 

Y n>N 

where K q = ^/(l/2 + J-)Cq and C q is as given by Corollary 1. 
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Proof. Recall Plancherel Formula E] |6fc(/)| 2 = \\f\\^. We consider two auxiliary functions g,h € 



feez 

H5, er ([— 1, 1]), given by the following formulae, 



g(x)= J2 hk W-ff' h & = E ^Tyf' W 
|fc|<[#] |fc|>[^]+i 

so that f = g + h on [—1, 1]. Moreover, we have ||p||2 < H/H2 and 
IWIl = E M/)| 2 < E (l + (M 2 )- s (l + (M 2 ) s |6 fc (/)| 2 <c- 2s ||/-/ [cA] ||^. (78) 

l fc l>[f]+! I fel >[#]+! 

We then proceed as in the proof of the previous theorem. The main point is the computation of the 
PSWFs series of g. We use the following computation of coefficients of the exponential e tXx . 

a n {e lX ) = J ^ e iXx i> n>c (x)dx = M<#«,c (^j ■ (79) 
It follows, by linearity, that 

1 , , ^ . , . (h% 



a n (h) = -=ti n (c) E 6*(/)^n,c(— J. (80) 
l fc l<[f] 

Then by using Schwarz Inequality as well as the fact that all quantities kir/c ar contained in the 
interval [—1,1], we get 



\a n (h)\ 2 < ^^||Vn, C 



E i^(/)i 2 E 1 1 ( 81 ) 

< Q + f W(c)\ 2 U n ,o\\l\\f\\l = Q + A„(c)||V„,c||^||/||| (82) 

One obtains 

1/2 

= I E i a «(5)i 2 ' 

\ n>N 



^ J( 1 / 2 + ^)T,\\^c\\Mc)\\fh, 



n>N 



which allows to conclude for (|75[) . Finally, by combining the previous inequality and the result of 
Corollary 1, we conclude for ([76]) . 

□ 

Remark 6. The last theorem may also be used for trigonometric polynomials bte lk7rx as long as 
their degree is smaller than c/tt, so that the second term in (|75p vanishes. 

Remark 7. We also have a bound of the error for ordinary polynomials. Indeed, if we consider the 
polynomial f(x) :— x J , then 

a n (f) = J y J ^ n ,c(y) dy = (-i) 3 'c~ 3 (j, n (c)i(i$ c (0), with i 2 = -1. 
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We can then use Proposition 5 to conclude that ifc 2 /\N < 1, then 

\\f-s N f\\l<c*Y.(^r] \^\ 2 - (83) 

This gives an alternative way to find and improve the results of Theorems 3 and 4, given in Jfij/. 
Indeed, the authors approach f by a polynomial, then consider coefficients of the PSWFs expansion 
for a polynomial. 

Remark 8. It is important to note that when the bandwidth c is comparable with the truncation 
order N, then the error bounds given by our previous two theorems are of the expected order 0(N~ S ), 
for the values of N in the decay region of the X n (c) . 



5 Numerical results 

In the first part of this section, we give a brief description of the numerical methods of computation 
of the PSWFs as well as the PSWFs series expansion coefficients of a function from the Sobolev 
space H s ([— 1, 1])- Note that Flammer's method is among the first methods that have been developed 
for the computation of the PSWFs inside [—1,1]. This method is well described in [5]. The explicit 
analytic extension of the PSWFs to the whole real line is given by D. Slepian. 

Recently, there is an extensive amount of work devoted to new highly accurate computational 
methods of the PSWFs, see [3] [TUJ [Til [22] • In particular, the methods given in [3j [26] are based 
on an efficient quadrature method on the unit circle that provides highly accurate values of the 
PSWFs inside [— 1, 1], as well as accurate approximations of the different eigenvalues £i n (c), n > 0. 
The methods developed in [TUl E] for computing the values of the ip niC (x) inside [—1, 1] and the 
eigenvalues /i„(c) are either based on an appropriate matrix representation of the finite Fourier 
transform operator Q Cl given by (|6"3"]) or a Gaussian type quadrature formula applied to Q c . 

In the sequel, we will adopt Flammer's method for computing the PSWFs. This choice is based 
on the facts that the values of the analytic extensions of the PSWFS over K, are easily obtained by 
this method. This later is briefly described as follows. We write as before the Legendre expansion 
of the PSWFs, 

OO ' 

VVcW = E PkM*)- (84) 

fe=>0 

oo ' 

where the sign means that the sum is over even or odd integers depending on whether the 

fc=0,l 

order n is even or odd. It is well known that the different expansion coefficients as we U as the 

corresponding eigenvalues X n (c) are obtained by solving the following eigensystem 

(fc + l)(fc + 2) „ 2an , , u/u , ^ 2fc(fc + l)-l , 2an 



A" +2 + (k( k + 1) + ";'_j c^i 



(2fc + 3)^(2fc + 5)(2fc + l) " +z K (2fe + 3)(2fc-l)' 

=c 2 /3[L 2 =X„(c)^, fc>0. 



fc(fc-l) JJ on 



(2k ~ l)^(2fc + l)(2fc-3) 

In [22], the author has shown that the analytic extension of the PSWFs outside [—1,1] is simply 
given by the following formula, 



<M») - EM)^ VfcTTT^ Jk+1 ^ cx) , vn>i, (85) 

IMnlCJI k>Q V CX 
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where 

, , 2n \E' k >o.i i k Vk + T/2 fit J k+1/2 (c) 

A*n(C) = — ; 

is the exact value of the n— th eigenvalue of the finite Fourier transform operator Q c 



(86) 



Remark 9. If f € Hp er , s > 0, then its different PSWFs series expansion coefficients can be easily 
approximated by the use of formula J79| ). More precisely for a positive integer K, an approximation 
a n (/) t° a n{f) is given by the following formula 

a*(f) = ^ E UMn,c (^) = a n (f) + e K , (87) 
where the bk(f) are the Fourier coefficients of f and where ej( = —^= l^n(c)bk(f)tpn,c I — I ■ 

V 2 |fc|>A'+l ^ ° ' 

Moreover, from the well known asymptotic behavior of the tp n ,c{x), for large values of x, see for 

example JlOf . one can easily check that ex = o ( j= — ^ — yi+i 1 ■ This computational method of the 

a n (f) has the advantage to work for small as well as large values of the smoothness coefficient s > 0. 

Remark 10. If f e H s ([-l,l]), where s > 1/2 + 2m, m > 1, is an integer, then f £ C 2m ([-1, 1]). 
Moreover since ip n . c € C°°(M), then the classical Gaussian quadrature method, see for example fSjj 
gives us the following approximate value &„(/) of the (n + 1)— th expansion coefficient a n {f) =< 
f,^n,c >, 

m 

5n(/) = *^2uJlf(xi)lp n ,c( x l) = a n{f) + Cn, (88) 
1=1 

with \e n \ < sup "' c . Here, b m is the highest coefficient of P m , and the different 

ve[-i,i] h m ( 2m r- 
weights uji and nodes xi, are easily computed by the special method given in J^j. 

Next, to illustrate the quality of approximation by the PSWFs, as well as to explain the contri- 
bution of the bandwidth c > in this quality of approximation, we give the following examples. 

Example 1: In this example, we show that the PSWFs outperforms the Legendre polynomials 
in the approximation of a class of functions from the Sobolev space H s ([— 1, 1]), having significant 
large coefficients at some high frequency components. To fix the idea, let A > 0, be a relatively large 
positive real number and let f\(x) = e lXx , x £ [—1, 1]. The Legendre series expansion coefficients of 
fx are given by 

o»(0)= / e lXx P n (x)dx = (i) n J — J„+i/ 2 (A). 



In this case, we have 



N 9 

|/a-5>«(°) P «II2 = - (^+i/ 2 (A)) 2 . (89) 

n=0 n>JV+l 



If c > is a positive real number, then the corresponding PSWFs series expansion coefficients of fx 
are simply given as follows, 

u n {c) = / e tXx ip n;C (x)dx = Hn(c)4> n , c {\/c). 
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In this case, the L 2 — approximation error is given by 



• / X 
c 



E N (c) = \\f-J2un(c)i>n,c\\i= £ \^{c)? (^c [ -) ) ■ (90) 

n=0 n>iV+l 



In the special case where c = A, the previous error bound becomes En (A) = |/i„(c)| 2 (i/j n _ c (l)) 2 . 

n>N+l 

Since VVc(l) = 0(y/n), and if we accept the numerical evidence that the exponential decay of the 
sequence (|/in(c)| 2 )„>o starts at [ec/4], then from ([89]) and ff90]) . one concludes that the PSWFs are 
better adapted for the approximation of the fx by its N— th order truncated PSWFs series expansion 
with c = A and N — [A]. More generally, if < c < A, then £ > 1 and the blowup of the tp n . c (£) 
implies that a n (c) = /j, n (c)ipn,c(^/ c ) has as lower decay that a n (X) = /in(A)Vvi,c(l)- Moreover, if 
c > A, then the decay of the |/x„(c)| 2 and consequently, the fast decay of the a n (c) is possible only 



if n lies beyond a neig hborhood of f > ^. This means that c = A is the appropriate value of 
the bandwidth to be used to approximate the function f\[x) — e lXx by its first N—th truncated 
PSWFs series expansion, with N = [A]. This explains the numerical results given in [25] concerning 
the approximation of the test function u(x) — sin(207ra), where the author has checked numerically 
that c = 207T is the appropriate value of the bandwidth for approximating u(x) by the PSWFs ip n ,c 
with a given high precision and minimal number of the truncation order TV. As another example, we 
consider the value of A = 50, then we find that 

50 50 

||/a - J2 a »(°)^ll2 « 3.087858£ - 01, \\fx~J2 a ™( 50 )^,5o||2 « 1.356604B - 08. 

n=0 n=0 

Example 2: In this example, we consider the Weirstrass function 

-!<,<!. (91) 

k>0 

Note that W s £ H s ~ e ([— 1, 1]), Ve < s, s > 0. We have considered the value of c = 100, and computed 

3 

W s ^n, the TV— th terms truncated PSWFs series expansion of W s with different values of — < s < 2 
and different values of 20 < N < 100. Also, for each pair (s, N), we have computed the corresponding 

1/2 



50 

- Yl (W s .N(k/50)-W s (k/50)) 2 



Table 1 lists the 



approximate L 2 — error bound Epj(s) 

; )\ i 

fc=-50 

obtained values of En(s). Note that the numerical results given by Table 1, follow what has been 

predicted by the theoretical results of the previous section. In fact, the L 2 — errors \\W S — HnW s \\2 

cc 

is of order 0(N~ S ), whenever N > N c ~ — — 1. In the case, where c = 100, N c — 67. The graphs 
of W 3 /4(x) and W 3 /^ N (x), N = 90 are given by Figure 1. 



Example 3: In this example, we consider the Weirstrass function 

x cos(2 fc 7ra;) , , . . 

/(») =S \hs y -1<«<1, * = l-4. (92) 

fc>0 

It is clear that / £ Hp er ([— 1, 1]), with \\f\l 2 , — fifes' ^ so ' we nx ^ ne va hie of the badnwidth c to 

fe>0 

100. In this case, we have ||/|| 2 ~ 1.0805838, ||/ - / [c/7r ]|| ff i « 1.203854. Next, we have computed 
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Table 1: Values of E n (s) for various values of n and s. 





s = 0.75 


s = 1 


s = 1.25 


s = 1.5 


s = 1.75 


s = 2.0 


n 


E n (s) 


E n (s) 


E»{s) 


E n (s) 


E n (s) 


E n (s) 


20 


4.57329E-01 


4.66173E-01 


4.85990E-01 


5.05973E-01 


5.23232E-01 


5.37227E-01 


30 


3.15869E-01 


3.11677E-01 


3.28241E-01 


3.48562E-01 


3.67260E-01 


3.82963E-01 


40 


1.06843E-01 


1.52009E-01 


1.91237E-01 


2.20969E-01 


2.43432E-01 


2.60523E-01 


50 


4.09844E-02 


6.88472E-02 


1.01827E-01 


1.26518E-01 


1.44809E-01 


1.58520E-01 


GO 


3.30178E-02 


2.09084E-02 


3.25551E-02 


4.28999E-02 


5.06959E-02 


5.65531E-02 


70 


3.15097E-02 


8.82446E-03 


2.51157E-03 


7.35725E-04 


2.33066E-04 


1.04137E-04 


80 


3.01566E-02 


8.55598E-03 


2.40312E-03 


6.87458E-04 


1.98993E-04 


5.80481E-05 


90 


2.67972E-02 


7.64167E-03 


2.14661E-03 


6.15062E-04 


1.78461E-04 


5.22848E-05 


100 


2.39141E-02 


6.72825E-03 


1.82818E-03 


5.10057E-04 


1.45036E-04 


4.19238E-05 



Figure 1: (a) graph of W s / 4 (x), (b) graph of W 3 /i t ^(x),N = 90. 



the different approximate L 2 — error E 



N 



r>0 



1/2 



^ ]T (/(A/50) ~ S N f(k/50)) 2 . In figure 2, we 

fc=-50 

have plotted the graphs of the actual error En versus the error bound of Theorem 5. Note that once 
N becomes larger than the critical value for the decay of the A n (c), which in our case, is given by 
N c = [ec/4] = 67, the theoretical error bound, given by Theorem 5 becomes very close to the actual 
error. 

Example 4: In this example, we let s > be any positive real number and we consider the Brownian 
motion function B s (x) given by as follows. 



fc>i 



-1< x < 1. 



(93) 



Here, Xj. is a Gaussian random variable. It is well known that B s € H s {[— 1, 1]). For the special 
case s = 1, we consider the band- width c — 100, a truncation order N = 80 and compute -B^at the 
approximation of B\ by its N— th terms truncated PSWFs series expansion. The graphs of B\ and 
Bi t N are given by Figure 3. 

Remark 11. From the quality of approximation in the Sobolev spaces H s ([— 1,1]) given in this 
paper and in J4^ \25f . one concludes that for any value of the bandwidth c > 0, the approximation 
error \\f — Sjv/||2 has the asymptotic order 0(N~ S ). Nonetheless, for a given f e H"([— 1, 1]), s > 



25 



Figure 2: (a) Graphs of the theoretical and the actual error bound for example 2. (b) Graph of 
the decay of A„(c). 



Figure 3: (a) graph of B\(x), (b) graph of Bi,n(x),N = 80. 

which we may assume to have a unit L 2 — norm and for a given error tolerance e, the appropriate 
value of the bandwidth c > 0, corresponding to the minimum truncation order N, ensuring that 
11/ — ^iv/lb < e, depends on whether or not, f has some significant Fourier expansion coefficients, 
corresponding to large frequency components. In other words, the faster decay to zero of the Fourier 
coefficients of f, the smaller the value of the bandwidth should be and vice versa. 
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